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We present the results on two-particle Bose-Einstein correlations measured in proton-proton col- 
lisions at center of mass energies of -Js = 0.9 TeV, y/s = 2.76 TeV and y/s — 7 TeV registered 
by the ALICE experiment at the Large Hadron Collider. Detailed analysis reveals that the three 
dimensional experimental correlation functions do not have Gaussian shape in two dimensions. 
We found that they are better described by an exponential functional form in the outward and 
longitudinal directions, while the sideward remains a Gaussian. This is interpreted as a result of 
a significant contribution of strongly decaying resonances to the shape of the emission region. 
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1. Introduction 

The Large Hadron Collider (LHC) at CERN started operating in the end of the year 2009 and 
has been colliding protons at the center of mass energies of y/s = 0.9 TeV, y/s = 2.76 TeV and y/s = 
1 TeV. At the LHC, the ALICE experiment ^ focuses mainly on measuring heavy-ion collisions. 
However, minimum-bias proton-proton collisions which provide the heavy-ion 'baseline' are also 
measured. The proton-proton data has been used by the ALICE [j^, |3j] and CMS [Q] collaborations 
in the two-pion Bose-Einstein femtoscopic analyses. In the detailed studies performed by ALICE 
on 0.9 TeV and 7 TeV data, in three dimensions, we showed that the shape of the experimental 
correlation function is clearly not Gaussian [ph. We tested different functional forms in each out, 
side and long directions and proposed a better fitting formula according to the best fit. In this work 
we extend the analysis to y/s = 2.76 TeV data. 

2. Data analysis 

Data samples of approximately 4 million proton-proton collisions at center of mass energies of 
y^ = 0.9 TeV, 20 million at y/s = 2.76 TeV and 60 million at y/s = 7 TeV were analyzed. The sub- 
systems used for the analysis were the Inner Tracking System (ITS), the Time Projection Chamber 
(TPC) and the VZERO detectors. All the events were required to have a reconstructed interaction 
point (primary vertex) within 10 cm from the center of the ALICE detector in the beam direction. 
The ITS and TPC were used for tracking in the pseudorapidity range \t]\ < 1.0 while the VZERO 
was used for a minimum bias trigger and to reject beam-gas and beam-halo collisions. Identifica- 
tion of pions was based on the energy loss of the particle information from the TPC. The analysis 
was always performed on primary particles and their selection was based on the minimum distance 
between the track and the primary vertex (so-called Distance of Closest Approach or DCA). Tracks 
were required to have DCA not greater than 0.018 + 0.035/?^ in the transverse plane and 0.3 cm 
in the longitudinal direction. We also applied specific procedures to suppress undesired two-track 
effects, such as splitting (one track reconstructed as two) and merging (two tracks reconstructed as 
one). For more details about the event and particle selection criteria see [^J. 
The analysis was performed in 8 ranges of total measured charged-particle multiplicity of the event 
N c h and 6 ranges of pair transverse momentum kj = 2 ■ 

3. Results of the femtoscopy analysis 

The femtoscopic analysis was performed in order to obtain the sizes of the particle emitting 
region. The Spherical Harmonics decomposition of the three-dimensional correlation functions 
was applied. This technique allows to represent the three-dimensional object as an infinite set of 
one-dimensional spherical harmonics functions. The symmetries of the pair distribution make most 
of the components vanishing and the first three of the non-zero ones, Cq, C\ and C\, capture most 
important information about the correlation [||, ^. The first one is the angle-averaged component, 
the second measures the difference between out and side while the third is the difference between 
transverse and long. The results show that the correlation functions are practically independent on 
the collision energy but they depend on the event multiplicity and (more strongly) on pair trans- 
verse momentum. Increasing the collision energy by an order of magnitude has less impact on the 
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experimental correlation functions than changing the multiplicity by 50%. All these effects are 
presented in Fig. |l[ 
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Figure 1: Comparison of the correlation functions for K + K + for a) different collision energies, b) different 
multiplicity ranges, c) different pair transverse momentum ranges. 

Let us focus on the C{j component. The femtoscopic effect, coming from the symmetrization 
of the two-pion wave function, is visible as the increase of the correlation function below q = 
0.5 GeV/c. The fall visible in the lowest bins is due to the Coulomb repulsion, but it has little 
impact on the Bose-Einstein peak and, therefore, on the extracted radii. We also see a broad peak 
extending in q from up to, at least, 1 .0 GeV/c. We performed a very detailed analysis of this 
effect in [||]. Because such structure can be seen in Monte Carlo data, which does not include Bose- 
Einstein correlations, it has clearly non-femtoscopic origin. Therefore, we can try to parameterize it 
using Monte Carlo simulations and use this parameterization to fit the data. From various formulas 
that we had tried to apply we found that it is quite well described by a Gaussian form in Cq and C\. 
Therefore, the final functional form for the background treatment is: 

(q-B m ) 2s 



B(q) = A/,exp (-q z A^) +B h e\p 



2Bl 



(3cos 2 (0)-l). 



(3.1) 



The analysis of the minijet origin of the non-femtoscopic background using angular correlations is 
described in \^. 

4. Fitting of the correlation function 

4.1 Gaussian fits 

The femtoscopic correlation function is defined by the Koonin-Pratt formula []8]] : 

C(q) = J S(r)\V(r,q)\d 4 r. (4.1) 

Usually it is assumed that the emission function describes a static source and has an ellipsoid 
Gaussian profile in space: 



5(r) =5(r ,r,,r/)~exp 
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where /?,- are the sizes of the source and r = [r uti r sidei r long\ is the pair separation vector. Using 
formula ( fO| ) and such parametrization of the emission function, one gets the following form of the 
correlation function: 

C(q) = C(q ou t,q S ide,qion g ) = 1 + Aexp {-R out q out -R side <l S ide - R iong a ion g ) > ( 4 -3) 

where X is the strength of the correlation (fraction of correlated pairs for which both pions were 
correctly identified). In the data analysis we must also take into account the Coulomb interaction 
between identical pions and strong interactions in the final state. However, since pions are emitted 
from sources which are on the order of 2 — 3 fm in size, the strong contribution is small and can be 
neglected []|]. The influence of Coulomb interaction is approximated using the Bowler-Sinyukov 
formula which assumes that the Coulomb part can be factorized out from the wave function of the 
pair and integrated separately. The modified correlation function, including Coulomb interaction, 
is given by the equation: 

C/(q) = (l-X)+XK(q) ■ [1 + Aexp (-R 2 0Ut q 2 0Ut - R 2 ide q 2 ide -RL g lL g )}, (4-4) 

where K is the Coulomb like-sign pion pair wave function averaged over the Gaussian source with a 
radius of 1 fm. Because the Monte Carlo generators do not include Bose-Einstein correlations, they 
can be used to estimate the non-femtoscopic background of the correlation function. In order to do 
this we fit the Monte Carlo simulated data with equation ( |Q| ) described in the previous section. 
The final form of the fitting function is therefore: 

C(q) =C(q out ,q s ide,qiong) = NCf(q out ,q s i de ,qi on g)B(q out ,q s j de ,qi ong ), (4.5) 

where N is the overall normalization. The Gaussian fits to the experimental correlation functions 
are presented in the left panel of Fig. f| It clearly underestimates the height of the peak. However 
the width of the correlation function is reproduced. 

4.2 Non-Gaussian fits 

In the previous section we have shown that the Gaussian fits do not perfectly describe the 
correlation functions obtained from the experimental data. In order to improve the fit we assume 
that the source function factorizes into out, side and long directions: 

5(r) = S(r out )S(r side )S(r long ). (4.6) 

This leads to the factorization of the correlation function itself: 

C(q) = i + XC(q out )C(q side )C(qi ong ). (4.7) 

We can change the functional form of each of the components of the source (and therefore the 
correlation function) independently. Three different functional forms of the source function were 
analyzed - Gaussian, exponential and Lorentzian. They have the desired feature, that the inte- 



gration in formula (11) can be performed analytically and lead to the Gaussian, Lorentzian and 
exponential forms of the correlation function respectively. It is reasonable to fit the data with func- 
tional forms other than Gaussian, especially in out and long directions, because resonances decay 
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after random time and this process is determined by the exponential decay law, which transforms 
into an exponential shape in space via the pair velocity (which by definition exists in out and long 
directions, and vanishes in side direction). We performed a study of all 27 combinations of the 
fitting functions for all the multiplicity and kj ranges. We found that universally in out direction 
the correlation function was best described by an exponential form, corresponding to Lorentzian 
emission function, which agrees with model expectations. In contrast, the side direction is equally 
well described by a Gaussian or a Lorentzian: we chose the former because the Lorentzian corre- 
lation function would correspond to exponential pair emission function with a sharp peak at 0. 
The non-Gaussian fit is shown in the right panel of Fig. ^[ One can see that the exponential form 
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Figure 2: Gaussian (a) and non-Gaussian (b) fits to the correlation function for the y/s — 7 TeV collisions 
with 23 < N^ < 29 and pairs with 0.3 < k T < 0.4 GeV/c. 

in out and long reproduces the data much better, for both the peak height and its width. 

4.3 Radii dependencies 

The dependencies of the femtoscopic radii on event multiplicity and kj extracted from non- 
Gaussian fits are shown in Fig. |3[ We stress that the radii cannot be directly compared to the ones 
from the Gaussian fits. We see that Ri ong and R s ide always fall with kj. On the other hand, R out has a 
different behavior - it firstly rises for the lowest kj ranges and then goes down. The radii also scale 
linearly with multiplicity. This scaling is observed in every kj range. Ru, ng and R out radii always 
grow with multiplicity while R out radius grows for the lower kj range and falls for the higher kj 
range. Gaussian fits show similar behavior. 

5. Conclusions 



We performed a detailed study of the Bose-Einstein femtoscopic correlations in proton-proton 
collisions at center of mass energies of y/s = 0.9 TeV, \fs = 2.76 TeV, y/s = 7 TeV collected by 
the ALICE experiment. We found that the correlation functions are almost independent on colli- 
sion energy, however they depend on the event multiplicity and on the pair transverse momentum 
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Figure 3: Pair transverse momentum and multiplicity dependencies of the femtoscopic radii from the non- 
Gaussian fits to the correlation functions from 7 TeV data. 

kj. The experimental correlation functions are clearly not Gaussians. We introduced an assump- 
tion that the source function factorizes, which leads to the factorization of the correlation function. 
Each component of the correlation function can be described by either Gaussian, exponential or 
Lorentzian functional form in out , side and long directions. From the fitting of all possible com- 
binations in all multiplicity and kj ranges we found that the exponential form in out and long and 
Gaussian form in side fit data best. The extracted radii from both Gaussian and non-Gaussian fits 
show similar behavior. 
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